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Abstract. It is shown that a string carrying a field of harmonic type can have 
circular vorton states of a new "throbbing" kind, for which the worldsheet geometry is 
stationary but the internal structure undergoes periodic oscillation. 

1 Introduction 

The purpose of the present work is to demonstrate the use of the formalism develped in 
a preceding article [1] for the treatment of fields in curved target spaces, by applying it 
to simple but non-trivial examples of the important special case of harmonic and other 
simply harmonious fields [21 13] on a string worldsheet. 

The present investigation will be restricted to geometric configurations of the sim- 
plest non-trivial type, namely circularly symmetric string loops in a fiat background 
spacetime, for which the metric will be conveniently expressible in cylindrical coordi- 
nates as 

ds^ = -dt^ + gMcp'^ + dg'^ + dz^ , (1) 

so that the string worldsheet will be specifiable by an expression for the radius g as 
a function of the time t at a fixed value of the longitudinal coordinate z which can be 
taken without loss of generality to be the origin z = 0. 

A systematic investigation of the dynamics of such a worldsheet has already been 
carried out [1] for conducting string models of the simple type for which the current has 
only a single degree of freedom, in the sense that the target space of the scalar field on 
the string is just one-dimensional - and therefor trivially fiat - the outcome being that 
if its energy is not too high the string will oscillate about a "vorton" type equilibrium 
state. Such a vorton state will be generically stable with respect to perturbations of the 
purely axisymmetic kind to which the present analysis will be restricted, but it has been 
shown [5] that they will be commonly, though not generically, unstable with respect 
to non-axisymmetric modes. It is to be expected that qualitatively similar behaviour 
will occur for scalar field models with more degrees of freedom [6], [7], so long as all the 
currents are generated by commuting symmetries of a fiat target space. 



1 



The novelty in the present work will be to consider the a situation of a qualitatively 
different kind that can arise when the relevant target space is not flat. Of course a 
curved target space might have no symmetries at all, in which case the currents in 
question would not even be conserved. The present work will however be concerned 
with the opposite extreme, in which the target space is highly symmetric, so that 
there will be many conserved current combinations, but with generators that do not 
commute. Attention will be focussed here on the simplest non-trivial possibility of this 
kind, namely the case in which the target space is just an ordinary 2-sphere, with metric 
ds^ = dX^ dX^ that will be expressible in terms of the usual coordinates X^ = 9 
and X^ = if hy 

dS2 = d^Vsin^^d^^ (2) 



2 Extrinsic motion of circular string worldsheet 

It will be convenient to describe the evolution of the worldsheet, within the background 
characterised by ([1]), in terms of unit timelike radial and spacelike transverse tangent 
vectors u'^ and u'^, and of a unit spacelike radial normal vector A^, that are given in 
terms of the coordinates x° = t, = 0, x'^ = g, x'^ = z by 

u^ = ^(6^ + g6^), u^ = -6^, X^ = ^(g6i; + 6^), (3) 

Q 

where a dot denotes differentiation with respect to the time coor dinate t and the 
Lorentz factor for the radial velocity g is defined as usual by 7 = 1/ a/1 — p^. 
The ensuing derivative formulae 

g g 

can be used to evaluate the second fundamental tensor as given [8] by the prescription 

K^S = r]:^^n.'. V^ = 77^^Va, (4) 
in which the first fundamental tensor of the worldsheet is specified as 

The second fundamental tensor of the time dependent circular worldsheet is thereby 
found to be ^ 

= l^^il^ Qu^u^ u^u^) . (6) 
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In the simple case for which the only external force is that of viscous drag by a 
static external background medium [H], which in this case will give a force density of 
the form 

F = f^' (7) 

with velocity dependent coefficient /, the corresponding equation of motion of the 
worldsheet will be given |H] in terms of the second fundamental tensor by an expression 
of the generic form 

T^'^K^J' = r , (8) 

in which T'^'^ is the relevant surface stress energy tensor, which will of course depend 
on the internal structure of the string. It can be seen that in this simple circular case, 
the ensuing differential equation for the radius g will take the form 

7^ T'^-'u^u, Q - ^ T'^'-u.u, = f. (9) 



3 Energy and angular momentum 

The invariance of the background ([1]) under the action of the time translation Killing 
vector k'^ and the rotation Killing vector g'^ defined by 

fc^ = = ^{u'- - , = 6^^ = gu^ , (10) 

( so that kyk^ = 1 and g^g'^ = g^) allows us to construct corresponding energy and 
angular momentum fiux vectors 

V>' = -k''%\ ji^ = Q'^Tf^ (11) 

which, subject to the variational field equations, will, as discussed in the preceeding 
work p, automatically satisfy the surface divergence conditions of the form 

V.p- = -F/^, v.j- = ^Mj^. (12) 

For an external force density of the postulated form (IT]) one thus obtains the work rate 
formula 

V,P" = 7^/, (13) 
and the angular momentum conservation condition 

V.J" = . (14) 
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It will be useful for what follows to rewrite these conditions in terms of internal 
worldsheet coordinates a*, with respect to which they will be expressible as 

V,r = 7^/, r = -7n^T/, (15) 

and 

%r = Q, j' = Qu^Tj\ (16) 
while the corresponding expression for the extrinsic equation of motion ([9]) will be 

-i^T^UiUjQ -^TiuiU, = f . (17) 
Q 

More particularly, with respect to the internal coordinate system that is induced on 
the worldsheet by taking a° = t, = 0, the corresponding expression for the intinsic 
metric of the worldsheet will take the form 

ds' = -i-d^' + ^)M(/.^ (18) 

and the corresponding expressions for the orthonormal frame vectors will be 

u'=^6l u' = -6l. (19) 
Q 

It can be seen that, with respect to these particular coordinates, the energy and angular 
momentum flux vectors will be given by 

Vi = -^^Tj J^ = Tj, (20) 
while the extrinsic equation of motion (ITTl) will be expressible more explicitly as 

^^T,,g-^J,, = f. (21) 



The total work rate formula ( USl) will take the form 

{QlV), = -QQf. (22) 
and the condition of the angular momentum conservation will take the form 

{-T:), = 0. (23) 

7 

It is to be remarked that fl2T]) can be used to eliminate the force magnitude / from 
fl22l) to give an intrinsic energy creation law of the form 

{gV), = gT,^. (24) 
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4 Generic harmonious case 



The formulae of the two preceeding sections are apphcable to classical string models of 
any kind . We now restrict attention to the harmonious characterised [1] by 

a Lagrangian L that depends only on the target space metric g^^ and the symmetric 
target space tensor defined - in the absence of gauge coupling, as will be assumed here 
- just by 

tt,--=r^X-X-. (25) 
This means that it its generic variation will have the form 

5L = ^^5X0^^ + ^^5g^ , (26) 

AB 

in which, as a Noether identity, we must have 

dL . dL 



dvoBc 



so that the coefficients will be specifiable by the expressions 

^ = ->^B ' ^ = ^""^ = ■ (28) 

AB 

in terms of the same symmetric target space tensor n^g . This tensor can be used to 
express the generic variation of the Lagrangian in the concise form 

6L = -'-n/SrVg\ (29) 

and to express the ensuing surface stress energy tensor as 

T^j=n,,X-X-^+Lg,^. (30) 



With respect to the coordinates of (ITTI) . using a prime for differentation with respect 
to 0, and a dot (as before) for differentiation with respect to t, we shall obtain 

th^B = J_ x-^'X^' - j^X^X^ , (31) 



and the stress energy components in (1201) will thus be given by 

1 ^ 

7' 



T,, = K,gX-X- - - L , T,, = n,gX-'X-' + L . (32) 
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Our investigation will be concerned with solutions that are axisymmetric in the 
strict sense [T] , meaning that the gradient fields X * and X-'^' are independent of (j), but 
not in the strong sense which would require that even the undifferentiated fields 
should be independent of (p. This means that X^' is allowed to be non-zero, but that 
we require X^' = X^" = 0. Under these conditions the total work rate formula (|22|) 
will take the form 

ig^^K^,X^X^-g^Ly = ggf, (33) 
and the angular momentum conservation law fl2Tl) will take the form 

{g^K,,X-'X-y = 0, (34) 

while the intrinsic energy creation law will be expressible in the form 

{gV^^AB^'' X^)- + ^K/ {g'ro,-')- = 0. (35) 

If the target space is only two-dimensional, and in particular if it is a 2-sphere as 
in the example dealy with in detail below, the complete system of dynamical evolu- 
tion equations will be provided just by the pair of internal equations flMl) and fl5^ in 
conjunction with the extrinsic evolution equation obtained by substitution from f l5^ in 
(I2TI) . Further input from the set of current pseudo-conservation laws constituting the 
complete system of internal field equations [H |2] will however be needed if the target 
space dimension is three or more. 

5 Quadratically and simply harmonious models 

Within the extensive category of harmonious models to which the foregoing formulae 
are applicable, a noteworthy subcategory is that of models that are quadratically har- 
monious, in the sense of being governed by a Lagrangian whose dependence on Ktg-'^ is 
just quadratic, so that it will be expressible in terms of fixed parameters m, n^,, 
in the form 

L = -m- ^K^Yo - ja^Yo'^ + j/3^Ki^''Kig^ , (36) 
with the usual notation rb = tb^-^, which gives 

An important special case is that for which = /J^,, so that the quadratic part is 
interpretable as a current cross product: this gives what is known as a baby Skyrme 
model [To], [11] when the target space is a 2-sphere, and it gives a fully fiedged Skyrme 
model [m [13] when the target space is a 3-sphere. 
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The quadratic special case for which = belongs to another noteworthy sub- 
category, namely that of simply harmonious models [2], which are characterised by a 
Lagrangian L that depends only on the scalar rb, as given by the formula 

m = , (38) 

so long as gauge coupling is absent, as is supposed here, so that with respect to the 
coordinates of f llip it will take the form 



'^ = ^^..^"'^^'-7^^..^^^^ (39) 



In this simply harmonious case we shall have 

I^AB = 1^9 AB , (40) 

with the coefficient k given by 

K = -2^. (41) 
dro 

In terms of this quantity, the intrinsic energy creation law f l35|) will be expressible in 
the form 

{qV^^' 9,,X^ ■ + K\gX-'X-'y = . (42) 



6 Minimally non-Abelian — spherical target — case 

Let us now concentrate on the minimally non-Abelian case, meaning that with the 
simplest non-flat target space geometry, namely that of a 2-sphere as given by ([2]). 
In such a case there are only two internal degrees of freedom,namely those of the 
independent field variables 6 and (p, their evolution will be fully determined just by 
the two preceeding conditions (l34l) and ( l35l) if the loop radius g is given in advance, as 
for example in the artificial case in which / is adjusted to hold the radius at a fixed 
value with ^ = 0. These two intrinsic evolution equations will also be sufficient, in 
conjunction with the total work rate equation ( !33|) if the external force magnitude / is 
given in advance, and thus in particular in the case of most obvious natural interest, 
namely that in which it is taken to vanish, 

f = 0. (43) 



The possibility of having configurations that, with respect to the rotation Killing 
vector g^, are symmetric not in the strong sense, which would require X-^' = but 
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in the less restrictive weak, albeit strict sense [T], as postulated here, depends on the 
existence of a corresponding symmetry in the target space, with generator V-* such that 

X^' = . (44) 

In the spherical case under consideration here, such a vector field could be chosen 
in many ways as a combination of the set of not just one but three independent target 
space Killing vector fields, for which the standard basis aa-^ is given [1] for a =1,2,3 
(corresponding to what are respectively interpretable as rotations about the West, East, 
and North poles) for = 9, = (p by 

= —sin ifi S-^ — cot 9 cos (p 5-^ , = cos (p 5-:^ — cot 9 sin >p 5^ , = 6^ . (45) 

(On planet Earth, in the roughly Jerusalem centered system favoured by cartographers 
since the time of Dante, the West pole is in the South Atlantic where the Greenwich 
meridian itersects the equator in the vicinity of the Gulf of Guinea, and the East 
pole is in the Indian Ocean, again on the equator but 90 degrees further East in the 
vicinity of the Bay of Bengal, while the North pole is of course in the middle of the 
Arctic Ocean. The places referred to in the Old Testament, including Jerusalem and 
particularly Noah's legendary landing place. Mount Ararat, are near the centroid of 
these three poles, opposite to what Dante called the antepode, which is about as far as 
possible from any major land mass, in the middle of the South Pacific.) 

There will be no loss of generality in choosing the coordinate system in such a way 
as to align V-^ with the last of these, that is to say with the generator of rotations about 
the North pole, which means that we shall have 

V^=nd^, (46) 

with a proportionality constant n that is evidently interpretable as a winding number, 
so that it must be an integer (which would have to be zero in the special case of strong 
symmetry). This simply means that the space gradients involved in the dynamical 
equations above will be given just by 

^' = 0, ifi' = n. (47) 

As explained in the preceding work [1], the existence of the target space Killing 
vector fields ( H5l) allows the two independent internal field equations to be expressed as 
conservation laws for the three currents given, for a =1,2,3, by 

Jai = nao^^g^^X^,, (48) 
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of which only two are independent. It is to be recalled that k is specified by the 
equation of state as a function of the quantity ro, which will be given in this case by 

rv = '^sin'e--f''0^ + sm^e^^). (49) 

It can be seen that the equation of conservation for the third of these currents, namely 

J,, = Ksm^e{^61 + n6l), (50) 

will take the form 

(^7Ksin2^(^)- = , (51) 

which contains just the same information as the equation ( l34l) for conservation of 
angular momentum, except in the special case of strong symmetry, n = 0, for which 
the angular momentum simply vanishes. It can also be seen that the only independent 
information obtainable from the conservation of the other two currents i and i is 
that of the internal energy creation equation (142|) . which will take the form 

[q'^^k^O^ + sm^e(p'')) ■ + K^n\siY?Oy = . (52) 

To obtain the complete system of equations of motion for the three independent vari- 
ables 6*, (p, and q, the internal dynamical equtions ( 1511) and (l52ll need to be supplemented 
by the information about the extrinsic motion that is contained in (l33ll which, in the 
force free case characterised by fH3l) . will take the form of the total energy conservation 
condition 

(^7^k(^2 + sin2^(^^) - ^7L)" = 0. (53) 



7 Harmonic separability for spherical target case 

It is evident that the preceeding set of three dynamical equations will immediately pro- 
vide two constant first integrals of the motion, namely a current conservation constant 

Q-i Ksm^e (p = <L , (54) 

obtained from (15T1) and a total energy conservation constant 

Q 73k,(^2 ^ gjj^2^ (p^)- Q-^L = <t, (55) 

obtained from (153|) . However the third dynamical equation ( 152|) will not be so con- 
veniently integrable in the generic harmonious case, for which the coefficient k is a 
variable function of the quantity tt) given by (H^ . 
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In order to proceed, we now restrict attention to the special case of a model that 
is not just harmonious but but actually harmonic, so that the coefficient k is just a 
constant. The harmonic case is characterised in terms of a pair of constants m and 
by a Lagrangian of the merely linear form 

L = —m^ — ^ K^rv , (56) 

which simply gives 

K= K^. (57) 

The presence of the (Kibble type) mass term is irrelevant for the purely harmonic 
equations (15T|) and (152|) that govern the internal fields on the world sheet if the latter 
is prescribed in advance, but for the actual calculation, via ( 153|) . of the evolution of 
the worldsheet the specification of m is indispensible, as it fixes the value of the string 
tension in the zero current hmit. 

In this special harmonic case, the first constant of the motion fl54p can be used to 
eliminate the variable (p, which will be given simply by 

c € 

g 7 sm 6 i^* 

and it is apparent that the second internal dynamical equation ( l52l) will also provide 
a constant ffist integral, which can be specified as the necessarily positive quantity 
given by the formula 

g'^m^ = , (59) 

using the notation tt)^ for the quantity obtained by changing the sign of the second 
term of the definition (149!) of tt), namely 

= ^sin2^ + 72(^2 ^sin2^^)2N /go) 
^2 



In the special harmonic case ( 1561) this notation can be used to rewrite ( 1321) as 



Too = ^{^^ + k'^* ' ^11 = (-^' + k'^* ' (61) 

and to rewrite the formula (l55i) for the energy constant in the form 

g'j{m'^ + ^K^tv^) = (E. (62) 

This leads to the discovery of a remarkably convenient separability property, whereby 
the external dynamical variable g can be decoupled from the internal field variables 6 
and 1^ by the elimination of iti^ between (!59l) and (134|) . The ensuing separated equation, 
for the radial variable g by itsel, can be seen to take the form 

SVW^ = m^^+^-. (63) 
z g 
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8 Throbbing vorton states 



It can be seen that the radial evolution equation (!63|) will give rise to an evolution that 
will be qualitatively similar to what has been found |4] for circular strings with just 
a single independent current variable, which is that the loop will oscillate periodically 
between finite minimum and maximum values of its radius g. 

More particularly, when the energy constant € is taken to have the minimum value 
compatible with a given value of the other constant, a^, a vorton type equilibrium state, 
with fixed radius 

Q=b, ^ = 0, (64) 

will be obtained. By minimising the left hand side of ( j63l) . it can be seen that such a 
vorton will be characterised by 

b2 = K^aV2m^ <E=^/2Klma. (65) 

Unlike the stationary, meaning strictly time independent (though not static, mean- 
ing strongly time independent) vorton states of the kind that familiar in cases when 
there is only a single current - or even when there are several currents if their generators 
commute - a vorton state of the non-Abelian kind considered here has the remarkable 
feature describable (using a term borrowed from the medical context of blood circula- 
tion) as throbbing. What this means is that, although the stress tensor and worldsheet 
geometry of the string loop are time independent, its internal fields undergo non sta- 
tionary oscillations. By substituting from ( l58l) and ( 159|) in ( 160|) . the field 9 can be seen 
to have a non-trivial time evolution given by 

b^^^^a^-cVsin^^-n^sin^^. (66) 

It can be seen that, whenever > -|- n^, the colatitudinal field 6 will oscillate 
symmetrically between a minimum, where 2v? ?,iv?6 = — va^-4r?c^, in the northern 
hemisphere, Q < 6 < n and a maximum with the same value of s\.v?9 in the southern 
hemisphere, n /2 < 9 < n. If > n^, such oscillating non-Abelian configurations can 
be viewed as perturbations of a strictly stationary single current configuration with 
longitude variable = n0 ± ct/b, for fixed equatorial colatitude, cos 6* = 0. For small 
amplitudes, the perturbations will have cos^^ oc cos{a;t} with = (c^ — n^)/b^. 

There can be no solution at all with < 2|nc|, but solutions of a rather weird 
kind will be possible for the intermediate range 2|nc| < < -|- n^, provided the 
supplementary condition < v? is also satisfied. For this parameter range, the field 
9 will be asymmetrically confined to a single one of the target space hemispheres, 
oscillating between values where 2n sin 9 = a^± Va^ - 4n2c2 without ever crossing the 
equator where 9 = 7r/2. Such oscillating non-Abelian configurations can be considered 
as perturbations of a strictly stationary single current configuration having longitude 
variable of the chiral form (p = n(0 ± t/b), with fixed colatitude 9 = arcsinA/|c/n|. In 
the small amplitude limit, the perturbations will have u"^ = 4(n2 — IncD/b^. 
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9 Harmonic separability for axisymmetric target 



The main motive for the preceeding work was to exhibit the behaviour of currents 
generated by target space symmetries that do not commute - and so cannot be made 
simultaneously manifest - by considering the simplest case for which non-commuting 
symmetries are present, namely that for which the target space is spherical. However 
it has turned out that whereas the one-parameter Abehan subgroup corresponding to 
axisymmetry is essential, the existence of the other non-commuting symmetries has 
played no qualitatively important role in the foregoing results, for which the indispensi- 
ble postulate was that the field equations in question should be not just harmonious but 
of the strictly harmonic form f lSB]) . which is all that is needed to obtain the constants 
and € as given by (j59p and ([22D in terms of the quantity 

ro' = ^g^^X-'X-' + ^'g^J-X-. (67) 

Elimination of this quantity then gives a separated radial evolution equation for g of 
exactly the same form (163) as for the special case of a target space that is spherical. 

For complete the separability of the system when the target space is two dimensional, 
it is sufficient that its metric should have the general axisymmetric form 

dS^ = q'^dw^ + d(^2 , (68) 

with q given as an arbitrary function of to. By setting uj = sin 9 it can be seen that 
this metric will take the spherical form ([2]) in the special case for which q = l/y/l — to^, 
and it will simply be flat if g = 1. The choice of the function q has no effect on the 
condition (|58|1 . which will simply go over to the form 

^ = — (69) 

and implementation, as before, of the strict but weak axisymmetry postulate, to the 
effect that to' = but (p' = n, for some non-vanishing integral value of n, will reduce 
(j67j) to the form 



to^ = -,uj' + j'e^' + ^ . (70) 

It follows that for all such cases there will be geometrically stationary throbbing vorton 
states characterised via ( |59|) by the same equations (!64l) and (!65|) as before, and thus 
with internal structure governed by a dynamical equation of the form 

f2-2-2 2 ^ 2-2 f^^^ \ 

b q zu = a — — — n to , (71) 



TO 



which is soluble by quadrature to give 



I bq^^__ (72) 
va2TO2 — — n^TO^ 
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The simplest example is of course the one provided by the case ^ = 1, namely the 
model having just a single complex scalar field, with amplitude w and phase (p, for 
which the target space is fiat, 

dS^ = dx^ ^ + dx^^ , x^=wcos(p^ X^=rosin(^. (73) 

This is the case for which the internal field model is purely linear, so that it will admit 
multiply conducting vorton states of the ordinary strictly stationary kind, in which the 
conserved currents are generated by the Abelian algebra of the target space translation 
group. This Abelian algebra is however just a subalgebra of the complete symmetry 
group: although the target space is fiat, its symmetry group is non-Abelian because it 
also includes rotations, which of course do not commute with translations. The presence 
of the conserved currents generated by such non-commuting rotations is what allows this 
familiar simple model to provide vortons not just of the usual strictly stationary kind, 
but also of the throbbing kind considered here. For the linear field model characterised 
by (173!) the corresponding quadrature (1721) with g = 1 can be evaluated exphcitly: the 
internal field amplitude w will throb in a manner given by the formula 



{o? + e'^cos{ujt}^ /2n^ , = Va^ - 4:cH'^ , u = 2n/b . (74) 



The concomitant formula for the throbbing of the internal phase variable (p can seen 
from (169|) to take the form 

ip = n(f) + arctan I ^"^^ tan| | i • (75) 



a2 + £2---i.2 

It follows that the Cartesian field components (1731) will be given by the expressions 

V 2v? = Va2+ ^2 cos{n0} cos{nbt} =f \/^ — sin{n0} sin{nbt} , 

\/ 2v? = Va^+ e"^ sin{n0} cos{nbt} ± Va^ — cos{n0} sinjnbt} , (76) 
(where the sign ± is that of the product cn) or equivalently by the complex combination 

+ = 77- exp{m0} (V 2nc exp{inbt} + — 2nc exp{—inbt}] . (77) 

When e"^ <^ (near the limits 2nc — )■ ia^) such a throbbing solution can be regarded 
as a perturbation of an ordinary stationary vorton configuration of the special chiral 
type, as given by (p = n(0 ± bt), with w = ^/\c/n\. 
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